Dynamical formation of a nonequilibrium subsystem under severe action 
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Formation of the nonequilibrium subsystem in dynamical processes during defect generation is 
simulated by means of molecular dynamics. A particular process of dissipation of the low-frequency 
acoustic emission into high-frequency equilibrium vibrations of lattice is studied numerically. Clear 
heuristic reasons are used to introduce different temperatures and entropies for equilibrium and 
nonequilibrium subsystems. Simple relaxation equation is proposed to describe time behavior of the 
nonequilibrium entropy. 
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I. INTRODUCTION 

Well known problem of the statistical physics and 
thermodynamics [l|, [2| is related to reducing of the re- 
versible laws of the micromechanics into irreversible laws 
of nonequilibrium thermodynamics 0, 0, Q ■ It invokes 
some simplifying assumptions which solve different par- 
ticular problems, but leads to a loss of generality. Af- 
ter many years of the study it still remains reasonable 
to seek for the new ideas and, maybe, return to initial 
sources of the problem. In solids, for example, with their 
complex internal structure the problem is related to ex- 
istence of weakly interacting subsystems. In many cases 
one can define relatively small spatially separated subsys- 
tems where an equilibrium state establishes quicker than 
in the whole system (so called, local thermodynamic equi- 
librium) [g, [3|. In another cases one can define spatially 
overlapped subsystems, which have their own tempera- 
ture and other thermodynamic parameters. For example, 
it takes place for the electron and phonon subsystems, 
which can be almost totally separated in the adiabatic 
approach [8[. The same can happen for the phonon sub- 
system and radiation Q, or in particular for acoustic 
emission [T(|- I n all these cases every subsystem may 
be in the equilibrium state, but complete system is not 
equilibrated. 

Formation of the defects in the solids under external 
impact provokes nonequilibrium processes. It transfers 
a part of the internal energy into thermal form. How- 
ever, this part does not transform into heat immediately 
and radiates in the first stages in a form of low-frequency 
acoustic waves. These waves gradually dissipate and sup- 
ply the energy to equilibrium thermal motion. As result, 
one can distinguish two weakly interacting subsystems of 
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high-frequency phonons (vibration modes) and the low- 
frequency acoustic emission. First subsystem is normally 
equilibrium before the external impact, and remains al- 
most equilibrium after it. The second one is created by 
the impact and strongly nonequilibrium. Present article 
is devoted to study of the interaction between these two 
equilibrium and nonequilibrium subsystems by means 
molecular dynamics simulation. 



II. MD-SIMULATION 

Let us regard to a particular configuration of the hard 
indenter pressing into 3D atomic specimen (Fig. [1]). The 
indcntcr (three upper atoms in the planar projection) 
moves with constant velocity 5m/s. The atoms possessed 
by lowest array are motionless. Behind the surface shown 
in the front of the picture 3D system contains 10 closed 
packed planes of atoms placed into HCP structure. All 
lateral faces are free. The sample consists from copper 
atoms interacting by Lennard- Jones potential. At some 
time moment the dislocations emerge in the atomic spec- 
imen. Fig. [T] presents exactly the stage when the dislo- 
cations start to enter into the system. It manifests itself 
by a bend of the curve for time depending total internal 
energy (curve 1, Fig. ^ or in the jumps of total po- 
tential energy (curve 2). Second subplot presents time 
dependence of the velocity for arbitrarily chosen particle 
of the system (which actually records its thermal mo- 
tion in curve 3). The thermal motion can be treated as 
high-frequency signal which is close to a harmonic one for 
separate interatomic (bond) vibrations. Its spectrum is 
severely oscillated in consequence of interference of large 
set of high-frequency phonons (region 2, Fig. (3^)- In 
the low-frequency region a large peak stands up sharply 
against the rest spectrum (peak 1). Performing low- 
frequency nitration over 100 time steps one can separate 
the low-frequency vibration modes (curves 4 and 5 in the 
last subplot). These modes gradually dissipate with time. 
Comparing curves 1-5 one can see directly that initiation 
of the low-frequency vibrational excitations coincides in 
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of the high-frequency vibrations, and seems to be enough 
to smooth of random fluctuations of the thermal motion. 
One can notice from the figure that equilibrium tempera- 
ture increases as energy of nonequilibrium low-frequency 
vibration transfers into equilibrium high-frequency mo- 
tion. As low-frequency wave or vibration motion differs 
from high-frequency one by time scale only, one can use 
this observation to define a nonequilibrium temperature. 
For a large enough numerical system (and certainly for 
the real solids) it is possible to take so long enough time 
average to suppress completely random fluctuations be- 
longing to the short range scales and separate large-scale 
perturbations caused by the entrance of the dislocations. 
Thus, digital filtration allows us separate equilibrium and 
nonequilibrium components from the initial data. 



FIG. 1: Geometry of 3D computer experiment positions of 
atoms at the 2900th time step. Time step is equal to 10.6 Fs. 



time with a generation of the dislocations and can be 
associated in this sense with strongly nonlinear acoustic 
emission. 

The external work performed during indentation in- 
creases total internal and potential energies (beginning 
parts of curves 1 and 2). Further, after t = 0.075ns, the 
potential energy decreases sharply (curve 2). A part of 
this energy transforms into defect energy and results in 
dislocation nucleation. Remaining part dissipates into 
heat motion via radiation in a form of the low-frequency 
vibration waves (curves 4 and 5 respectively). These 
acoustic waves scatter on high-frequency phonons and 
as result are damped exponentially. This process of the 
nonequilibrium state relaxation gradually leads to slow 
increase of equilibrium temperature and entropy. 

During indentation the entrance of the dislocations re- 
peats from time to time (see points t = 0.125ns; 0.175ns; 
0.26ns in Fig. Due to this some fraction of the low- 
frequency nonequilibrium wave packets always presents 
in integral thermal motion. Low-frequency component 
of wave motion differs from high-frequency phonon field 
by spatial and time scales only. So, to describe it one 
can use the same language as for the ordinary phonons. 
In particular, thermal motion of the phonons is charac- 
terized by temperature and entropy. The temperature 
can be treated in both manners: as mean energy per one 
degree of freedom (averaged over ensemble) or as mean 
energy of the vibrations of each particle (averaged over 
time) [ill, To get equilibrium velocity of the par- 

ticles one can subtract a record of nonequilibrium low- 
frequency vibrations (curve 4, Fig. [2]) from the complete 
one (curve 3). 

To control the procedure let us recalculate remaining 
energy the spectrum. It is shown in Fig. [3)d), where 
the peak 1 corresponding low-frequency vibration packets 
vanishes. Such a curve in Fig. 2] is plotted after average 
over interval of 100 time steps contains about 20 periods 



III. PHENOMENOLOGICAL MODEL 

Let us now apply phenomenological approach to de- 
fine nonequilibrium entropy. Kinetic equation for pro- 
duction and annihilation entropy in nonequilibrium sub- 
system may be written as follows [HI, [l4| : 



at = ls{T ° 



T st - Tis), 



(1) 



Here To is temperature of internal sources producing 
nonequilibrium heat at, so 7s To is an intensity of the 
entropy sources for the nonequilibrium subsystem, T st is 
temperature of internal acceptors of the dissipation for 
nonequilibrium heat and "fsT st is power of entropy sinks 
for the nonequilibrium subsystem, T±s is relaxation term. 
The value T s t has stationary meaning of nonequilibrium 
temperature to which the system tends during evolution. 
The last term in EqQ] describes transfer of the entropy 
from nonequilibrium subsystem to equilibrium one due 
scattering of low-frequency vibrations on high-frequency 
phonons. 

The validity of the above assumptions has been 
checked using another example of acoustical emission at 
formation of corrosion microcracks into bulk of welded 
seam. As a result of numerical simulations the distribu- 
tion of number of acoustic emission signals n on ampli- 
tude A during a kinetic stage is found to be described by 
Gaussian distribution [151 ]: 



n{A) = ngexp 



(A -(A)) 2 
2a 2 



(2) 



where no is a normalization, (...) is averaging operator, o 
is a standard deviation. The square of signal amplitude 
is energetic parameter per one act of acoustic emission. 
So, the relation (2) coincides with standard Gibbs dis- 
tribution which standard deviation has a meaning of the 
temperature T — 2a 2 /ks- In a stationary state one can 
assume T — T st , where T st is stationary temperature in- 
volved into kinetic equation JTJ). 
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FIG. 2: Time scanning: 1 total (internal) energy; 2 total potential energy; 3 velocity of arbitrary particle; 4, 5 low-frequency 
record of particle velocity after filtration. 




f, GHz f, GHz 

FIG. 3: Spectrum of total record (a) and difference one (b): 1- law frequency peak; 2 - main region of spectrum 



The same relations have been recently restored from 
another example of acoustic emission appearing during 
high-temperature deformation of aluminum. The acous- 
tical emission is caused here by Porteven-le Shatelie effect 
related to a short-run coherent slip of large family of dis- 
locations which form single slip trace and to a spatial 
correlation between deformation processes in mesoscopic 
scale [ill [TtJ • Temperature dependence of activation vol- 
ume the has been derived in the Ref. [H[ in the form: 

V = 2/o + Aexp (— J , (3) 



where y = 0.47±0, 12, A = 0.0027±0, 0004, t = 130±39. 
If the temperature is high enough one can neglect addi- 
tive term yo, expand the exponential over temperature 
T r and reduce the dependence ([3]) to the effective Gibbs 
distribution: 

n = n o ex p(~]^-f) > ( 4 ) 

Here n is number of separated pulses of acoustical emis- 
sion, forming given region of activation volume, e is ac- 
tivation energy of acoustic emission and the constants in 
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Thus, constant t has a meaning of inverse activation en- 
ergy. 



t, ns 



FIG. 4: Time dependence: 1 - equilibrium temperature; 2 
nonequilibrium temperature 



the equations and ^ are related as follows: 



I = /;,,< rp j —j^jT ) , 4=-^-. ("») 



To summary: using molecular dynamic simulation we 
have investigated a mechanism of nonequilibrium subsys- 
tem formation due to dislocation defect generation. This 
process produces low-frequency nonequilibrium phonons 
(acoustic emission) which dissipate into high-frequency 
equilibrium lattice vibrations. The temperatures of equi- 
librium and nonequilibrium subsystems (as well as conju- 
gated to them equilibrium and nonequilibrium entropies) 
have been introduced on a base of clear heuristic rea- 
sons and illustrated using some additional examples of 
the acoustic emission. Simple relaxation equation is pro- 
posed to describe the nonequilibrium entropy evolution. 
A comparison of the results with experimental studies 
of acoustical emission on another structural scale is per- 
formed. 
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